Quantum Phase Tomography of a Strongly Driven Qubit 
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The interference between repeated Landau- Zener transitions in a qubit swept through an avoided 
level crossing results in Stiickelberg oscillations in qubit magnetization. The resulting oscillatory 
patterns are a hallmark of the coherent strongly- driven regime in qubits, quantum dots and other 
two-level systems. The two-dimensional Fourier transforms of these patterns are found to exhibit 
a family of one-dimensional curves in Fourier space, in agreement with recent observations in a 
superconducting qubit. We interpret these images in terms of time evolution of the quantum phase 
of qubit state and show that they can be used to probe dephasing mechanisms in the qubit. 
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Superconducting Josephson devices recently emerged 
as a platform for exploring coherent quantum dynamics 
in solid state systems [1]. Due to their macroscopic di- 
mensions, these devices feature strong coupling to RF 
fields [1, S S] 5 t)e used to study new quantum 

phenomena associated with strong driving such as Rabi 
oscillations in the multiphoton regime [Sl, lSK Landau- 
Zener- Stiickelberg- type (LZS) interference 0,BQ, Bloch 
oscillations [l^, and qubit-photon dressed states 

In the LZS regime 0, H, 0, the qubit undergoes re- 
peated Landau-Zener (LZ) transitions at a level crossing, 
with adiabatic evolution between crossings [12]. Inter- 
ference between subsequent LZ transitions results in an 
oscillatory dependence of qubit magnetization in the fi- 
nal state on the detuning from the level crossing and the 
driving amplitude 0, [l^. The LZS effect is related to 
earlier observations of photon-assisted transport in quan- 
tum dots 3, m and in superconducting systems [16, 17], 
which exhibit multiphoton sidebands with oscillatory de- 
pendence on RF field amplitude. Although the observed 
oscillations washed out more quickly at high RF power in 
those devices than in the qubits 0, H, @] , in all cases the 
oscillations originated from the LZS interference effect. 

A new regime of strong driving was reported in a re- 
cent work [l8[, in which a qubit was driven through a 
manifold of several states spanning a wide energy range 
up to 120 GHz. The observed LZS interference indicated 
that even such strong driving was not detrimental for 
coherence. The change in qubit magnetization induced 
by the driving pulse exhibited complex checkerboard-like 
patterns in the two-dimensional phase space parameter- 
ized by the DC magnetic flux and RF driving ampli- 
tude. These patterns displayed a multiscale character, 
with multiphoton resonance lines on the finest scale and 
LZS interference fringes on a larger scale, and additional 
complexity due to resonance and interference effects in- 



volving several pairs of energy levels during each pulse. 

In an attempt to better understand the observed 
patterns, the authors of Ref. [18] employed a two- 
dimensional Fourier transform (FT). Unexpectedly, 
the FT revealed a highly ordered structure of one- 
dimensional arcs joined together to form lemon-shaped 
ovals in Fourier space, in contrast to the familiar Bragg 
peaks in the Fourier images of periodic patterns. Most 
surprisingly, these arcs were found to connect the high 
and low wavenumber regions, which are associated with 
the multiphoton resonances and LZS interference fringes. 

In this article we explain the lemon-shaped structures 
observed in [l^], first using a quasiclassical phase argu- 
ment and then fully microscopically. We analyze the FT 
of the transition rate (see Figjl]), which can be measured 
using short excitation pulses [9]. Then we consider the 
FT of qubit magnetization produced by long pulses. 

Our analysis reveals a relation between the lemon- 
shaped structures and the coherent dynamics of the 
qubit. In fact, because the Fourier transform inverts the 
energy variable and maps it onto the time variable (scaled 
by h), we find that the lemon arcs can be interpreted as 
an image of the time dependence of the quantum phase 
of the qubit (see Eq.([2|) below). This relation, as we shall 
discuss, can be exploited to probe fundamental aspects 
of qubit dynamics such as decoherence and dephasing, 
and to measure the decoherence times T2 and T2*. 

The intensity of each point on the curve in Fourier 
space arises from a particular ray-like section of the LZS 
pattern (see FigHj), with the section direction in one-to- 
one correspondence with the time interval beween level 
crossings. The section- by-section mapping to Fourier 
space is reminiscent of tomographic imaging, and realizes 
a "tomogram" of the time evolution of the qubit phase. 

Employing the FT to image quantum phase is famil- 
iar from the work on the mesoscopic Aharonov-Bohm ef- 
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FIG. 1: Tomographic imaging of qubit phase evolution. The 
pattern of LZS oscillations in (a) the transition rate (|12|) and 
(b) its Fourier transform, which exhibits a family of parabolic 
arcs dSI) forming lemon-shaped ovals. Along each arc, ke and 
kA represent the time separation t2 — ti and phase gain 0^2, 
Eq.(|5|), between subsequent level crossings, respectively. The 
Fourier intensity at each such point is mapped from the re- 
gion near a ray s/A = in the {A,e) plane, where the pa- 
rameters A, s all yield the same time interval and phase gain 
between level crossings. Four of these rays and the corre- 
sponding points in Fourier space are shown in matching col- 
ors. A Sawtooth-like driving signal (inset) was used, with the 
decoherence rate r2 = ^a; in ([T2)). 



feet 20|, 21 1, whieh used the dependence of conductance 
on magnetic field. In our approach, however, the time 
dependence of the phase is reconstructed using a two- 
dimensional FT where the axis associated with the de- 
tuning from qubit level crossing plays the role of time. 
We also note that in recent work 22, 23, \2M a to- 



mographic reconstruction of the Wigner function on the 
Bloch sphere was performed. The procedure used in 
Refs. 22, [2^, 24 1, which employs controlled rotations of 



the qubit state following its Rabi oscillations, is differ- 
ent from that used in the present work. Our image of 
qubit time evolution is obtained in Fourier space. Also, 
because of the nature of the LZS effect, it only provides 



information about the relative phase of the qubit |0) and 
|1) states, rather than the entire Bloch vector. 

To emphasize aspects common to different experiments 
that have used harmonic 0,(1, 0, [ll|, [13] , sawtooth- 
like [13], and bi-harmonic 0] RF driving, we consider 
a generic periodic driving of the qubit: 



h (h{t) 



Mt). (1) 



where h{t) is the energy detuning from an avoided cross- 
ing, periodically modulated by the driving field g{t) = 
g{t -\-T) with amplitude A and zero mean, f g(t)dt = 0. 
For simplicity here we focus on the case when g{t) has 
one maximum and one minimum per period. 

Away from the level crossing, the qubit evolves adia- 
batically as a superposition of the states |0) and The 
LZS interference can be expressed through the relative 
phase of the states |0) and |1) gained between subsequent 
passages through a level crossing: 

(t){A,e) = j h{t)dt = e{t2 - h) - A J g{t)dt. (2) 
ti ti 

The times ti^2 of level crossing are the solutions to 

Ag{t) = £ {gmin < e/A < ^max), (3) 

represented graphically in the inset of Figjl^ by the inter- 
sections between lines of fixed detuning and the driving 
signal. Quasiclassically, the LZS contrast can be modeled 
by a sum of functions cos(^(A, e)), one per each choice of 
ti,2 in ([2|). We consider a position-dependent wavevector 



{kA, k,) = ±{VaHA V,(/)(A, £)), 



(4) 



where ± accounts for the contributions of e^^^^^^'^^ . Eval- 
uating the derivatives in (|4]) and, noting that the net con- 
tributions of \/sti^2 and Va^i,2 vanish due to (|3]), gives 



{kA,h) = ±{-(j){2,t2 - ti), (/)• 
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g{t)dt. (5) 



Crucially, Eq.([5|) defines a curve parameterized by a sin- 
gle variable u = e/A^ which is the only parameter upon 
which the times ti^2 found in Eq.(|3]) depend. 

We illustrate this mapping by an example of sawtooth 
driving ( Fig [T] inset) with g{t) linear between the points 



5(nT) = -5((n±i)r) = l. 



(6) 



In FiglH the points in Fourier space arising from differ- 
ent s/A sections are denoted by dots with colors match- 
ing those of the corresponding rays £/A = u'm the (A, e) 
plane and of the constant detuning lines in the inset. The 
ks and kA coordinates of those points correspond to the 
time separation (Fig. la inset) and the phase gain, Eq.©, 
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FIG. 2: Transition rate {[12]) and its FT for sawtooth- like 
driving for r2 = four times larger than in Fig[T] A 3D 
projection plot of Fourier intensity is shown to illustrate the 
exponential decay WrrikA^ks) oc e"^^^^. 

between subsequent level crossings. In this way, the 
curves (j5j) reproduce the time evolution of qubit phase. 

As shown in FiglU each ray maps to a family of points 
([5]). The reason for this multiplicity is that, besides the 
sign ± in Eq. ([5|), Eq.(|3]) has multiple solutions t'l = 
ti + niT and ^2=^2 + for each e and A, where T 
is the period of driving and ni,2 are arbitrary integers. 
Because § g{t)dt = 0, all such solutions yield the same 
phase gain same value of kA- However, the 

corresponding values of k^ are displaced by (77,2 — ni)T, 
generating the periodic family of arcs shown in Fig{T)3. 

Another class of solutions to Eq.(|3]) describes subse- 
quent passages through the level crossing in the same 
direction: t2 = ti -\- nT. In this case ti is unconstrained 
and, because zero phase is gained over a single driving pe- 
riod, we obtain a discrete set of points (/c^, ks) = (0, nT) 
irrespective of e, A. As displayed most clearly in Fig|2t, 
the FT intensity indeed exhibits peaks at ks = nT. The 
peak positions k^ = nT agree with the spacing huj be- 
tween multiphoton resonances in the (A, s) plane. 

To find the form of the curves in Fig [1)3, we solve Eq.© 
for the case of sawtooth driving, Eq.®. Without loss of 
generality we select -T/2 < ti < < ^2 < T/2 and find 

t2 = -ti = r/2, r = T{A- e)/2A. (7) 

Evaluating the phase ^f2 = It^ 9{t)dt = i(l " £^/A^)T, 
we obtain parabolic arcs in Fourier space: 

{kA, ks) = ± (-(1 - r/T)r, r + nT) , (8) 

< r < T, where the term nT was added to k^ to account 
for the multiple solutions to Eq.([3]) discussed above. 
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FIG. 3: Test of the result ()15|) for the Fourier intensity dis- 
tribution in the lemon, a) FT of the transition rate in Fig|2] 
after the image was doubled by extending W to negative A 
as W{—A^£) — W{A^e). b) Intensity in a mid-lemon section 
(blue line) compared with a similar section through a lemon 
in FT of the undoubled image in Fig|2](red line). 

Similarly, in the case of harmonic driving, the solutions 
of Eq.([3]) are ^2 = —ti = ^ diTccos{s / A) ^ which gives the 
phase ^f2 = cos{ujt)dt = {2/u)y^l - e'^ /A'^ . Substi- 
tuting these results into Eq.© we obtain the sinusoids 

ujkA/2 = ±sm{ujks/2) (9) 

which were observed in Ref. [18]. 

Now we turn to a microscopic analysis of qubit dy- 
namics based on the Hamiltonian ([T]) to which we add 
classical noise to model decoherence: h{t) = h(t) + fe(t). 
The transitions between qubit states |0) and |1) can be 
analyzed most easily in a rotating frame where 

« = -f(A.% T)' ^W=Ae-*«, (10) 

with ^{t) = Jq h{t')dt' . Perturbation theory in A yields 
the rate of transitions between the states |0) and 

A 2 r rt+5t 

W= lim — / / {e-''^^^'^e''^^'^^)ssdtidt2, (11) 

5tT2'>l A5t J 

where r2 = ^ is the decoherence rate. We average over 
5e{t) using the white noise model: l^e^^^^'^'^)-^^^^'^^)) = 
e-r2|ti-^2|^ where d(l){t) = j'^ 6e{t')dt' . 

To find the rate W in closed form, we use the Fourier 
series e^^^'^^ = e*^^ /me~*^^S where the coefficients 
fm can be expressed through the error function of com- 
plex argument for the case of sawtooth driving, or Bessel 
functions for the case of harmonic driving Q . Using the 
appropriate Fourier series in (pTj) and performing the in- 
tegration over ti and ^2, we obtain the expression 

W(s A) = — V (12) 

m= — 00 ^ / 
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At cj ^ 27rr2 this expression describes non-overlapping 
resonances (see Fig(T]), while at uo ^27rr2 it describes the 
partially dephased regime of Ref. H (see FigE]). 

To evaluate the Fourier transform of the transition rate 
WpTikA, ks) = JJZo e-'^^^-''^^W{e, A)dedA, it is con- 
venient to return to expression (pT]) . Because the phase 
(j){t) = et — Ag(t')dt' is linear in e as well as in A, we 
can easily bring the Fourier transform of (pT]) to the form 

r rt+5t 

aJJ^ S{k, + ti - t2)S{kA + (}){2y-^^^''-'^Utidt2 

with a = /S?{2tiY /4.6t and (t){^ defined in (p. This result 
can be simplified by performing the integration over t2 
with the help of the delta function 5{k^ + ti — ^2), giving 

WFT{kA,ke) = |A2a;e-r^l'=^l /(5(fc^ + 0?2)dti, (13) 

where t2 = ti -\- ks. Result ([13]) illustrates the effect of 
dephasing on the lemon structure through the prefactor 
e^^^l^el (^gee Figl2t), which arises from the exponential 
decay in time e~l^2-^il/^2^ ^^^^^ consistent with the in- 
terpretation of ks as a time variable. 

It is instructive to compare this behavior with the ef- 
fect of ensemble averaging, modeled by random offsets Se 
with a gaussian distribution. Because the phase factors 
in (pT]) are linear in £, the ensemble- averaged FT is 

(H^fT(fcA,fce))e„sOce-r^l'==le-^^'=S A=(fe2). (14) 

Through this dependence, instrinsic dephasing and en- 
semble averaging, i.e. T2 and T2*, can be distinguished. 

The lemon boundary obtained from (p!3|) for a generic 
g{t) agrees with the quasiclassical result (|5|). Indeed, the 
range of kA for which FT intensity is nonzero, at a fixed 
ks^ are determined by the extrema of the function 
ti. Writing 5t^(l)\2 = 9(^2) — 9 {h) = we recover Eq.(|3]). 

For the case of harmonic driving, g{t) = cosuot, we 
can evaluate (p!3|) by noting that = {sm{u;t2) — 
sin(a;ti))/a; = {2/ij) sin{^ijks) cos{uj{ti + ^k^)). The in- 
tegral over ti in (p!3|) then yields 

WFTikA^ ks) = , (15) 

2^^^\n\\u;ks)-k\ 

for < ^1 sm.{^uoks)\^ and zero elsewhere. We see that 
WFT{kA^ks) is concentrated inside the region bounded 
by the sinusoids ([9]) with square root singularities at the 
boundary. Similar behavior with a square root singular- 
ity in FT intensity is obtained for the sawtooth case, as 
illustrated in FigO Because Eq. (p!5|) is derived with the 
FT taken over — 00 < A < 00, the LZS pattern in Fig|2] 
had to be doubled to obtain the FT in Fig|3l 

Finally, lemon structures are also exhibited by the FT 
of the qubit steady-state population. The lemon arcs 
with multiple periods, clearly visible in FigHJ arise be- 
cause of a nonlinear dependence of saturated population 



a) b) 
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FIG. 4: Qubit magnetization and its FT. Shown is the magne- 
tization of a qubit driven to saturation: m — (Fi — F^)/ (2VK + 
Ti + r;), r; = Tie-^/'^BT^ ^YiQiQ Vi (r;) is the down (up) 
relaxation rate [9^. Results are shown for sawtooth driving 
with parameter values: decoherence rate V2 — tempera- 
ture k^T = 1.2/icj, relaxation rate Fi = 8 • 10~^a;, frequency 
zy = 400 MHz, level splitting A = 12 MHz. 

on W, with quadratic nonlinearity giving double period, 
cubic nonlinearity giving triple period, etc. This multi- 
plicity of periods was also observed in the data [18]. 

In conclusion, FT-based tomography of two- 
dimensional LZS patterns is a general technique 
that offers a way to image the quantum phase evolution 
of qubit s and other quantum systems. In the simplest 
case of a driving signal with just one maximum and 
one minimum per period, we predict a chain-like lemon 
structure in Fourier space which is in perfect agreement 
with recent observations. 
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